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For systems of lattice anyons like Majorana and parafermions, the unconventional quantum statis¬ 
tics determines a set of global symmetries (e.g., fermion parity for Majoranas) admitting no relevant 
perturbations. Any operator that breaks these symmetries explicitly would violate locality if added 
to the the Hamiltonian. As a consequence, the associated quasi-degeneracy of topologically non¬ 
trivial phases is protected, at least partially, by locality via the symmetries singled out by quantum 
statistics. We show that it is possible to bypass this type of protection by way of specifically engi¬ 
neered gauge fields, in order to modify the topological structure of the edge of the system without 
destroying the topological order completely. To illustrate our ideas in a concrete setting, we fo¬ 
cus on the Ze parafermion chain. Starting in the topological phase of the chain (sixfold ground 
degeneracy), we show that a gauge held with restricted dynamics acts as a relevant perturbation, 
driving a transition to a phase with threefold degeneracy and Z 3 parafermion edge modes. The 
transition from the Z 3 to the topologically trivial phase occurs on a critical line in the three-state 
Potts universality class. Hence, to all effects and purposes, the gauged Ze chain realizes the Z 3 
parafermion chain. We also investigate numerically the emergence of Majorana edge modes when 
the Ze chain is coupled to a differently restricted gauge field. 

PACS numbers: 05.30.Rt, 11.15.Ha 


I. INTRODUCTION 

The pursuit of topological quantum information pro¬ 
cessing with mesoscopic arrays^^— of Majorana chains^ 
is grounded on well investigated solid state technologies 
and showcases a new bridge between condensed mat¬ 
ter physics and engineering. With current experimen¬ 
tal techniques, naturally available topological quantum 
matter is ill suited for quantum information processing)^ 
Hence, one should try to design synthetic topological 
quantum matter with enhanced desirable characteristics. 
Mesoscopic physics naturally sets the stage for this pro¬ 
posal because the microscopic understanding of interact¬ 
ing topological quantum matter is incomplete. Hence 
one would like to exploit naturally occurring topological 
states, not as they come but as building blocks in meso¬ 
scopic arrays. 

The Majorana chain itself for example is not natu¬ 
rally available, and so the proposals for its realizatiort^ 
are good examples of the quest for synthetic topologi¬ 
cal quantum matter. The parafermion chain^ is an even 
more ambitious proposal along these lines;^ motivated by 
a problem that parafermions cannot, at least directly;^ 
wholly resolve in the end: the limited computational 
power of Majoranas. With parafermions the ide a^^i^^ is 
to modify mesoscopically the edge of a fractional quan¬ 
tum spin Hall insulator by placing alternating s-wave 
superconducting and ferromagnetic islands in order to 
obtain, at low energies, an effective Id lattice system of 
hybridized parafermion zero-energy modes. This basic 
mesoscopic blueprint has excited theoretical interest in 
various types of parafermion chains, with varying degrees 
of exoticism in their quantum phase diagramaH^— and 


also varying potential for mesoscopic realization^^ 

Parafermions, and in particular Majoranas, are exam¬ 
ples of topological zero-modes. Such modes are good for 
quantum information processing for three basic reasons: 
they are local, isolated in energy, and stable. As a conse¬ 
quence, the manifold of ground states of the system that 
hosts these modes inherits desirable properties. First, its 
degeneracy is stable. Second, it is possible to perform ro¬ 
tations in the ground manifold by coupling local probes 
to the topological modes. And third, there is an energy 
gap to non-topological excitations that sets the time scale 
for these manipulations. 

In spite of their crucial practical value, the stability 
of topological zero-energy modes poses a problem for in¬ 
cremental design. Take for example the basic blueprint 
for parafermions of the previous paragraph. The quan¬ 
tum statistics of the emergent lattice parafermion field is 
characterized by an i nteg er p >2, and on the mesoscopic 
setups of Refs. [Ill and[ij for example, p can only take the 
values p = 2m, with m an odd, positive integer. As a con¬ 
sequence, topological qutrits (corresponding to p = 3) are 
not included in this otherwise very attractive proposal. 
This is a problem for the circuit model assisted by topo¬ 
logical information processing, since quantum algorithms 
with qutrits are the best studied ones after qubits. Now, 
since the ground degeneracy of the p = 6 parafermion 
chain is six, it is natural to seek for a design that mod¬ 
ifies the basic platform minimally, just enough to halve 
the degeneracy. But the degeneracy is protected by the 
stability of the topological zero-energy modes! It should 
be clear from this example how, in general, incremental 
design principles can be at odds with “topological” ones. 

In this paper, we will propose a method based on sym- 


2 


metry manipulation without symmetry breaking to mod¬ 
ify in a controlled fashion the structure of the set of 
zero-energy modes at the edge of topological quantum 
matter. Our method in a nutshell starts by identify¬ 
ing discrete gauge fields naturally present in mesoscopic 
setups and naturally coupled to protecting symmetries. 
After providing dynamics to the gauge field, some or all 
of the topological degeneracy disappears. However, if the 
gauge fluctuations are properly manipulated, it is possi¬ 
ble to remove only some of the topological degeneracy. 
As a consequence, the topological structure of the edge 
changes in a controlled and predictable way: some zero- 
modes disappear and the remaining ones control the new 
ground manifold of the system. 

Given the symmetry analysis underlying our idea, one 
may ask why not just break explicitly the symmetry in 
question. An explicit symmetry breaking term would 
in general obtain the same result as that obtained by 
the gauge field without breaking the symmetry, and, in 
one space dimension, these two scenarios are often con¬ 
nected by a duality transformation. The answer is rooted 
deeply in the structure of anyon fields. As we will show, 
due to anyon statistics, the symmetries that must be ma¬ 
nipulated in our framework cannot be explicitly broken 
without violating locality. Effectively, this observation 
provides a symmetry-based explanation of the stability, 
or “protection,” of the zero-energy modes. By coupling 
the symmetry to a gauge field rather than breaking it, we 
manage to modify topological zero modes without break¬ 
ing locality. This idea evolved from our investigation of 
the effect of phase slips on the Majorana chain,— since 
they obtain a Z2 gauge field minimally coupled to the 
Majorana degrees of freedom. 

The outline of the paper is as follows. In Sec. m we 
explain our general ideas using a familiar example, the 
Majorana chain of Kitaev. The key point is that a I12 
gauge field arise naturally in physical realizations of topo¬ 
logical chains in terms of charged fermions. Making this 
local Z2 gauge field dynamic by mesoscopic manipulation 
allows to modify the topological edge structure without 
explicitly breaking fermion-parity. In the following sec¬ 
tions, we generalize these ideas to the parafermion chain. 
In Sec. uni we recall the basic theory of parafermions, 
discuss the parafermionic generalization of the Majorana 
chain and show that the global symmetry protecting 
the parafermionic zero-modes cannot be broken explic¬ 
itly without violating locality. In Sec. IlYl we introduce 
the gauged parafermion chain where parafermions of de¬ 
gree p are naturally coupled to a Zp gauge field. By 
means of an exact duality transformation we show that 
by restricting the dynamics of the Zp gauge field to a 
nontrivial subgroup of Zp we can obtain parafermions 
of reduced degree without explicitly breaking the pro¬ 
tecting symmetry. We show evidence that our construc¬ 
tion works as predicted by investigating the period of 
the supercurrent in a gauged parafermion ring junction 
as a function of the strength of the gauge fluctuations. 
Restricting the gauge held dynamics to a Z2 subgroup. 


we obtain a supercurrent periodicity consistent with Z3 
parafermionic edge modes. In SectionIVl we conhrm in 
another way that we have Z3 parafermionic zero modes 
in this gauged Zg chain by investigating the phase dia¬ 
gram of the system with the DRMG algorithm. The criti¬ 
cal exponents of the transition between the topologically 
non-trivial and trivial phase correspond to the univer¬ 
sality class of the three-state Potts model for any hnite 
strength of restricted gauge fluctuations, that is, on a 
critical line. We also compute the phase diagram of the 
model with effective Z2 parafermionic zero modes (Ma- 
joranas). In both cases, the modified topological zero 
modes are stable and exist in an extended region of cou¬ 
pling space. We conclude with a summary and outlook 
in Sec. ED 


II. DESTROYING THE EDGE MODES OF THE 
KITAEV MODEL WITH GAUGE FIELDS 

In this section we briefly recall the gauged Majorana 
chai n^^i^^ in order to keep a concrete example in mind 
as we proceed to investigate the more complex case of 
parafermions. 

Let us introduce two Hermitian Majorana operators 
Qj, bj per site site j = 1,..., L of a chain, such that 

— 0, — ‘ 26 ij. 

Majoranas are related to Pauli matrices, 

i=i i=i 

and fermionic creation and annihilation operators, 

CLi — Ci -t- , (2) 

-i6j =-Ci-I-cj = Ci(-l)"'*-b cj. (3) 

Let F = denote the total fermion number oper¬ 

ator. Then the operator of fermionic parity 

L 

(- 1 )^ = ( 4 ) 

i=l 

measures the total number of fermions modulo two. 

The Majorana chain of Kitaev^ is specified by the 
Hamiltonian 

L L-l 

ILk — '\h ^ U ^ ^ bjdj-^-ij (5) 

i=i j=i 

describing the competition between a topologically non¬ 
trivial phase for \h/J\ < 1 and a trivial phase for the 
complementary regime. The points \h/J\ = I are critical. 
In its topologically non-trivial phase, the Majorana chain 
displays twofold quasi-degeneracy protected by fermionic 
parity (— 1 )^. This symmetry cannot be broken explicitly 
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due to a superselection rule enforcing vanishing matrix 
elements between sectors distinguished by fermionic par¬ 
ity. Such a superselection rule can be traced back to the 
requirement of locality of the physical Hamiltonian. 

It has been shown in Ref. El that the quasi-degeneracy 
of the topologically non-trivial phase cannot be lifted by 
local perturbations. This result suggests that there is no 
way of modifying the topological structure of the edge by 
local perturbations. But this inference is in fact incorrect 
for a physical realization, as opposed to the ideal repre¬ 
sentation, of the Majorana chain.— The reason is that 
the Majorana chain will typically be realized in terms of 
electrically charged fermions, e.g., electrons, which carry 
the electromagnetic 17(1) gauge structure. Induced su¬ 
perconductivity breaks most of that symmetry, but a 
Z 2 part survives— This residual gauge symmetry, corre¬ 
sponding to an invariance of energy levels under a local 
multiplication of the Majorana operators aj, bj by a fac¬ 
tor of —1, is no longer manifest in the Majorana chain 
due to a tacit but specific gauge choice. Restoring the 
gauge fields normally used to absorb the effects of these 
residual gauge transformations leads to the Hamiltonian 

L L-l 

Hgk = -\h'^atbi - ^ [\Jbjaj+iTj + ktj] . ( 6 ) 

Z=1 j = l 

of the gauged Majorana chai n^^i^^ The ^2 gauge field is 
represented by Pauli matrices , rj, placed on the links 
j = 1,..., L — 1 of the chain. For k ^ 0, the gauge field 
is quantum dynamic— 

The gauged Majorana chain must be supplemented 
with gauge constraints identifying physical states. Since 
the 17(1) symmetry originates from electromagnetism, all 
physical states should obey the Z 2 descendant of Gauss 
law. The chain has local symmetries 

Gi = (-iai6i)rf, Gl = Tl_,{-\aLbL), (7) 

G,=Tf_,{-\a,b,)Tf (j = 2 ,...,L-l), ( 8 ) 

and enforcing the Gauss law corresponds to demanding 
that the local symmetries Gj for j = 2,..., L of the 
gauged Majorana chain act as the identity on physical 
states, Gj = 1 — 

The perturbation proportional to h in Eq. dHl) plays a 
subtle but definite role in tying the fate of the topolog¬ 
ically non-trivial phase to the gauge field. The reason 
may be traced back to the effect of the gauge field on 
the topological boundary modes of the Majorana chain. 
For h = 0, the gauged Majorana chain displays exact 
zero-energy modes 

[77gkU=o, ai] = 0 = [Hgh\h=o,bL], (9) 

perfectly localized at the end points of the chain and com¬ 
pletely insensitive to the gauge field. This insensitivity 
disappears as soon as h ^ 0. In the absence of gauge fluc¬ 
tuations, K = 0 , it is still possible to compute explicitly 


the evolution of the edge modes with h. We obtain 

7/=ai+^(j) aiTf_i...Tf, ( 10 ) 

i^2 

L-l , 

7r=^L + Xl(7) (11) 

with 

[77gkU=o,7;] = 2 ij(j)'^ 6 LrE_i...Tf, ( 12 ) 

[77gkU=o,7r-] =airE...rE_i. (13) 

The gauge field and Majorana (quasi) zero-energy bound¬ 
ary modes are now explicitly intertwined. In general, any 
gauge invariant perturbation, no matter how small, will 
suffice to intertwine the topological zero modes with the 
gauge field. 

For some physical setups, due to the stiffness of the 
phase of the superconducting condensate, the dynamics 
of the gauge field is strongly suppressed, k « 0. In this 
case, the gauge fields may be removed from the theory 
by a specific gauge choice in order to recover the topolog¬ 
ically non-trivial Majorana chain ([S]). However, we have 
showed in a previous publicationii that a finite value 
for K may be obtained through quantum phase slips in 
a mesoscopic implementation of the Kitaev model. At 
finite k, the properties of the Kitaev model are funda¬ 
mentally altered, because k measures the strength of a 
relevant perturbation of fermionic parity. Any non-zero 
value of K will remove the topological degeneracy of the 
Majorana chain for large enough system size. But why 
is this the case? 

One explanation comes from solving the gauge con¬ 
straints d?]) for Tj, giving 

j . ^ 

Tj = Y[i-\akbk) = (-1)^ II (14) 

fc=l k=j+l 

In this way the dynamics of the gauge field rj drives non¬ 
local parity measurement of the part of the chain to the 
left and to the right of site j, which is, in general, incom¬ 
patible with the topological degeneracy. While the gauge 
fields are local degrees of freedom, they have highly non¬ 
local effects on the electronic degrees of freedom, which 
makes them prime candidates for the manipulation of 
topological phases without breaking protecting symme¬ 
tries. Moreover, the effects we are describing are not tied 
to the mean field approximation. Very recently there has 
been some progress in the investigation of particle num¬ 
ber conserving Majorana chains^SiSi and, as should be 
clear from our arguments, we would not expect our con¬ 
clusions to be any different for the Z 2 gauged versions of 
such models. 

To summarize, discrete gauge theories can emerge in a 
natural, controllable way in physical realizations of topo¬ 
logical zero-energy boundary modes. Moreover, for the 
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Majorana chain, one can destroy the topologically non¬ 
trivial phase phase by gauging the global symmetry of 
fermionic parity and making the corresponding Z 2 gauge 
field dynamic. For parafermions, a richer structure is 
possible. This is the subject of the rest of the paper. 

III. PROTECTION OF THE EDGE MODES OF 

THE PARAFERMION CHAIN BY LOCALITY 

In this section we will briefly recall the basics 
of parafermions in one dimension, introduce the 
parafermion chain as a natural generalization of the Ma¬ 
jorana chain, and explain how the ground degeneracy 
of the parafermion chain is protected by the interplay 
between locality and anyonic statistics. There is a very 
recent review on parafermions Ref. 0 that the reader may 
consult for further details on implementation proposals. 

Parafermions are described in terms of a natural gen¬ 
eralization of the Majorana operators of the previous sec¬ 
tion. A one-dimensional, lattice-regularized parafermion 
field, defined on sites i = has the following 

properties^ 


PiAj = ujAjTi {i < j), (15) 

TiTj^wT^ri ii<j), (16) 

AiAj = ujAjA^ (i < j), (17) 

Tf = 1 = (-l)P-iAf, (18) 


where uj = ^ p-th root of unity. The or¬ 

der p = 2,3,... of the parafermion field characterizes 
the exchange (Eqs. (I15L (TlBl) . and (HD) and exclusion 
(Eq. (HD) statistics of these degrees of freedom.— In ad¬ 
dition, the parafermion field is unitary, 

r,ri = i = A,A^ (19) 

Up to a phase, parafermions of order p = 2 are Majo¬ 
rana operators. More precisely, 

1^1 ~ (P ~ ^)- ( 26 ) 

It is in general true of everything that follows that setting 
p = 2 will recover some standard facts about fermions. 

Following the reasoning for the construction of the Ma¬ 
jorana chain by pairing the second parafermion on each 
site with strength J with the first one the next site, one 
obtains the parafermion chain^ 

E(rlA, + H.c.) - i ^\a,f1+i + H.c.){2l) 

^ i=l ^ i=l 

For h = 0, the parafermion chain displays exact zero- 
energy parafermion edge modes, 

[i7p,ri]=0=[77p,Ai], (22) 

perfectly localized on the end points of the chain. As 
a consequence, the ground energy level is p-fold degen¬ 
erate, and, in fact, every energy level is at least p-fold 


degenerate. This is the topological degeneracy of the 
parafermion chain. Small values of will not remove 
the topological (quasi-)degeneracy of the ground energy 
level, but may remove the topological degeneracy of ex¬ 
cited energy levels 4*^ 

For the Majorana chain there is a strong link between 
topological degeneracy and fermionic parity. Is there 
something similar for parafermions? In order to obtain 
a precise answer to this question it is necessary to relate 
the parafermion field to a charged field and an underlying 
Fock space. Let us introduce creation and annihilation 
operators Cj , Ci satisfying the relations^ 


C'|^’ = 0 = C'f, 

(23) 

c\cy= t, 

(24) 

CiCj = ujCjCi (*<j). 

(25) 

C,C] = QC}C, (z<j), 

(26) 


(1 = l,...,p—1 ) where the bar denotes complex conjuga¬ 
tion, u) = The number operator for these “Fock 

parafermions”— is 

2m—1 

(27) 

and the Fock vacuum satisfies (7^10) = 0, Vi. Let P = 
denote the total parafermion number operator. 

Then 

ujP = (28) 

measures the parafermionic parity. 

The analogy between and (—1)^ is compelling, 
but it is not clear yet whether is a symmetry of the 
parafermion chain. The relation between Fock and stan¬ 
dard parafermions is^S 

V,=C, + cl^^-^\ (29) 

A, + (30) 

It is illuminating to compare Eqs. @ and (USD, and Eqs. 
dSD and (1501) . Since then 

rj A, = + cr^cl (31) 

Hence, 

a;^ = nrlA,. (32) 

i=l 

It is straightforward now to check that the parafermionic 
parity is indeed a symmetry of the parafermion chain. 

Similar to the fermionic parity (—1)^ in the Majo¬ 
rana case, the parafermionic parity distinguishes the 
various ground states of the parafermion chain. The 
parafermion field transforms under this symmetry as 

= uiTi, AiUj^ = ujAi. 


(33) 
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The stability of the parafermionic zero modes is ex¬ 
plained by the fact that is always a symmetry of any 
local Hamiltonian of parafermions: for these systems, 
the global Zp symmetry cannot be explicitly broken 
without violating locality. 

Why is this so? In order to answer this question, it is 
useful to introduce an explicit realization of parafermions 
analogous to the relation between Majoranas and Pauli 
matrices, Eq. O- First, associate to each site of the 
chain a pair of unitary p x p matrices Ui,Vi, often called 
clock or circulant matrices, that commute on different 
sites and otherwise satisfy 

= ojU,V,, Vf = Ul = 1. (34) 

We can represent a single pair of clock matrices as 



/I 

0 

0 • 

■ 0 

\ 



1 

0 • 

■ 


0 

UJ 

0 • 

■ 0 


0 

0 

1 • 

• 0 


0 

0 

■ 

■ 0 



0 

0 

0 • 

• 0 

u = 






, ^ = 






0 

0 

0 • 

■ 0 



0 

0 

0 • 

• 1 


Vo 

0 

0 • 

• UjP- 

V 


Vi 

0 

0 • 

• 0/ 


Then, 

r^ = V,U^-l...Ul, A, = r,f/*. (36) 

In this representation we see explicitly the string respon¬ 
sible for the non-trivial exchange statistics. Now it is 
possible to show that operators Oi localized in a finite 
region around site i, meaning that they are constructible 
as a product of parafermions in that region, commute for 
disjoint regions if and only if they are invariant under the 
symmetry of modular conservation of parafermion 
number. As a consequence any local Hamiltonian of 
parafermions is invariant under this symmetry as well: 
this symmetry cannot be explicitly broken without break¬ 
ing locality. 

Let us point out in closing an interesting relation be¬ 
tween Fock parafermions and electron fractionalization. 
For p = 2m with m odd, the Fock parafermion field is 
literally^! the mth rood of the electron field. More ex¬ 
plicitly, the composite fields satisfy canonical 

anticommutation relations, 

{C'r,C™} = 0, = (37) 

Tunneling between electronic and parafermionic systems 
can then be modeled by terms of the form —T/lC™ -I- 
H.C., where / is a canonical fermionic annihilation op¬ 
erator. Practical applications of this observation were 
explored in Ref. [13. 

For the mesoscopic setups of Refs, [llj and [13, it is 
indeed the case that the integer m characterizes the frac¬ 
tionalization in the system, since the quasiparticles of 
the fractional spin Hall insulator in the mesoscopic ar¬ 
ray carry charge e/m. The superconductors on the edge 
induce pairing, and the relation p = 2m is then a reflec¬ 
tion of the fact that it takes 2m quasiparticles to form a 


Cooper pair. It is interesting to notice that the m = 1 is 
allowed. In this simplest case, there is no fractionaliza¬ 
tion. The 2d insulator is just a spin Hall insulator, and 
the parafermion field becomes (up to a phase) a lattice 
Majorana field. Ref. [3 offers an up to date overview of 
proposals for realizing parafermions in mesoscopic arrays. 

IV. CHANGING THE EDGE MODES OF THE 
PARAFERMION CHAIN WITH GAUGE FIELDS 

In this section we will show how to selectively remove 
some of the zero-energy edge modes of the parafermion 
chain. By analogy to the gauged Majorana chain, the 
basic idea is to first restore the gauge symmetry of the 
parafermion chain, and then allow for suitably engineered 
gauge fluctuations in order to split in energy some gauge 
sectors and not others. The sectors that remain unsplit 
support parafermionic topological edge modes of reduced 
order. 

The generalization of the Z 2 gauge field of the Majo¬ 
rana chain is the Zp gauge field r^, ai satisfying 

T'i = = 1 ) ( 38 ) 

and commuting otherwise. This algebra is the same alge¬ 
bra as for the clock matrices Ui,Vi of the previous section, 
but we use a different notation here for clarity. Then the 
Hamiltonian of the gauged parafermion chain is 

h ^ 

H,p = -^J2^rlX + H.c.) ( 39 ) 

L-1 

- E [ 2 + H.c.) + + H.c.) . 

i=l 

In general, the design of the mesoscopic setup meant to 
realize the parafermion chain provides control over the 
gauge fluctuations, characterized by k and an integer a = 
0,... ,p — 1. The values k = 0 and/or a = 0 correspond 
to quenched gauge fluctuations, that is, no fluctuations 
at all. For a > 1 we say that the gauge fluctuations are 
restricted. The reason will become clear as we proceed. 

The local Zp symmetries of the chain are 

Gi = (rlAi)al, Gl = (40) 

G, = a,.iiTlA,)al (z = 2,...,L - 1), (41) 

They seem a priori independent of the details of the gauge 
fluctuations, and, being local, they cannot be sponta¬ 
neously broken)^ What is then the meaning of the of 
the integer parameter a? The answer is that local gauge 
fluctuations characterized by a > 1 generate highly non¬ 
local measurements of parafermion number modulo p/a 
which split the ground state multiplet, leaving behind a 
topological degeneracy of reduced degree. This can be 
seen by removing explicitly the gauge redundancy and 
performing a duality transformation, as we will show 
next. Later in this section we will show that this re¬ 
duced ground state degeneracy corresponds to topologi¬ 
cal parafermionic edge modes of reduced degree. 
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A. Changing topological degeneracy with 
restricted gauge fluctuations 

Let us start by rewriting the gauged parafermion chain 
in terms of the clock matrices of Eq. (135]) , 

= + (42) 

2^1 

by virtue of the transformation Eq. O- The result¬ 
ing model is akin to, but not quite the same as a Zp 
Higgs model if a ^ 1. In terms of clock matrices, the 
gauge symmetries are Gi = Ui(j\, Gt = ai-iUiO-j for 
i = 2,..., L — 1 and Gn = ctn-iUn- The global symme¬ 
try uj^ simply reads 

L 

io^=l[U,. (43) 

2 = 1 


Eor K = 0, the gauge-hxed parafermion chain iL|p 
is just the celebrated clock model (see Ref. [2^ for an 
up-to-date review), and one could translate back to 
parafermions to recover the standard parafermion chain. 
But for K ^ 0, the gauge-fixed Hamiltonian shows ex¬ 
plicitly how local gauge fluctuations at finite k affect 
the parafermionic degrees of freedom in a highly non¬ 
local way. The strings in iJ|p are not local neither in 
terms of clock or parafermion degrees of freedom. On 
the parafermions, they correspond to nonlocal measure¬ 
ments of parafermion number modulo p/a. 

To see why this corresponds to a reduction of the 
ground state degeneracy, we note that the gauge-fixed 
Hamiltonian may be brought into a local form by means 
of the duality transformation 

2 

(i = l,...,L), 

(z = 1,...,L-1), (50) 

Ul ^ Vl 


By introducing simultaneous eigenstates |{uj}) of the 
clock matrices Uj, 

U,\{uj]) (j = l,...,L), (44) 

and simultaneous eigenstates |{crj}) of the clock matrices 

aj, 

(j = 1) (45) 


exact for finite chainsi^ The resulting Hamiltonian is 

L -1 

2 ^'^ ' ^ 2 


= - |(t2 + H.c.) - I + H.c.) (51) 


L-1 




j>i 


we can define gauge-fixed states on the tensor-product 
space of the matter and gauge fields. They are of the 
form 


|{wj})gf = !{%}) ® \{aj = ^(ai - u*)}), (46) 

i>3 

with integers ai defined for i = 2,...,L. The states 
belonging to the gauge sector labeled by the at obey the 
relation 


and the global symmetry is mapped onto a boundary 
symmetry, uj^ i —. 

The Hamiltonian H|p’^ makes the role of the gauge 
fluctuations particularly clear. Let us focus on the phys¬ 
ical sector {ai = 0) for simplicity. Then, if we ne¬ 
glect finite-size corrections, that is, the boundary term 
hVL/2 -I- H.C., H|p’^ takes the simpler form Hi, with 


H, = - 




-f H.c. (52) 


Gi\{uj})gf — uj *|{uj})gf (i — 2,..., L), (47) 

and the physical sector corresponds to the choice = 0. 
Projecting the Hamiltonian onto the gauge-fixed states 
leads to the gauge-fixed Hamiltonian 

4^|p = -^E(t4. + E/) (48) 

^ 2=1 

- E [f +1 ’ 

i=l j>i 

where the operator Ui acts in the natural way on the 
gauge-fixed states, 

( 49 ) 


Let a/p = a!/r, with a' and r relative primes, and sup¬ 
pose for simplicity that p = rs with r and s also relative 
primes. Then, the equivalent of Eq. (l43l) with the 
matrices Ui of the dual model ((CT|) , gives 

= E“, (53) 

since = 1. Hence, the “longitudinal held” k 

breaks the symmetry but not completely if a > 

1: The Zig symmetry 0 ;”^° survives the gauge fluctua¬ 
tions unscathed. Conversely, the symmetry is 

completely broken. In this way, gauge fluctuations ob¬ 
tain speciflc relevant perturbations that reduce but do 
not remove completely the topological degeneracy of the 
parafermion chain. 
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B. Effect of restricted gauge fluctuation on 
topological boundary modes 


Restricted gauge fluctuations obtain a bulk relevant 
perturbation that partially removes the topological de¬ 
generacy of the parafermion chain. What is the effect 
on the edge modes? Here we address this question in 
terms of the Josephson effect, generalizing the celebrated 
anomalous periodicity of the supercurrent in the Majo- 
rana chain for investigating gauged Parafermion chains. 

Let us take p > a = s > \ for concreteness, so that 
a/p = 1/r, with r and s relative primes. Then the exact 
zero modes 


fi = r^, Ai = 


(54) 


of the gauged Parafermion chain with h = 0 realize 
parafermions of reduced degree s. These topological 
modes commute with the symmetry that is effec¬ 
tively, but not explicitly, broken by the restricted gauge 
fluctuations. Hence, we expect these zero modes to sur¬ 
vive the gauge fluctuations because they live in sectors 
of constant and will not become split in energy. 

We can make this picture quantitative and very close 
to the situation for the gauged Majorana chain if we re¬ 
call a technical result about clock matrices. Let us write 
for the clock matrices of order q = r,s. Then 
we have the identifications^^ 


with n the unique solution o^ 


(55) 


(r -|- s)n = 1 (mod p). 


(56) 


In particular, 

cr® ^ (erf ^)® (g) (57) 

Since^ 

(rf f"(fTff® = e'2"/’'(Tf ®, (58) 

we can make the simpler identification 

(59) 


In terms of the identification of Eq. (15^ , we can rewrite 
the gauged Parafermion chain in reduced form. 


^gp = -^E(rlA. + H.c.) (60) 


The site-dependent numerical phases are s roots of 
unity obtained from diagonalizing the r) ^. Since, in par¬ 
ticular, af* = 1 , they can be removed by a gauge choice 
in the definition of the parafermion field. Hence the situ¬ 
ation is as follows: A parafermion chain coupled to a hp 
gauge field with restricted fluctuations {a = s = p/r) is 


equivalent to a parafermion chain coupled to a gauge 
field with unrestricted gauge fluctuations (a = 1). Based 
on our experience with the gauged Majorana chain, we 
expect that the topological boundary modes 

Ai = A^ Al = AI (61) 

will be gapped out by the Zr gauge field, leaving behind 
the Zs zero modes of Eq. (IMl) . 

We can check this claim explicitly by extending the 
theory of the 47 r Josephson effect in the Majorana chain 
to our gauged parafermion chain, or more precisely, to 
iLpg to avoid unnecessary numerical cost. To drive a su¬ 
percurrent across the gauged parafermion chain it is con¬ 
venient to close the system into a ring threaded by a flux 
$. The Hamiltonian for this parafermion ring junction is 

iJ'($) = - ^(d+ i7.c.) (62) 

(the Fock parafermions Ci, Cj of charge 2e/p were in¬ 
troduced in the previous section). The tunneling term is 
also coupled to the gauge field in order to preserve gauge 
invariance. The gauge symmetries of the ring junction 
are 

G'i = (T^Ai) (erf V (* = * + L)- (63) 

Let us consider for concreteness a Z 2 gauge field cou¬ 
pled to parafermions of order p = 6, so that r = 2 
and s = 3. Fig. [IJa) shows the six lowest energy lev¬ 
els of i7gp($) in the Z 2 gauge invariant sector and for 
K = 0 = h. Any gauge sector would show the same 
result since they are all degenerate in energy. The 
full 127r Josephson effect expected of the usual p = 6 
parafermion chai n^^i^^ is clearly visible. For k/J = 0.2 
and h = 0, the period of the supercurrent has decreased 
to 67 r, c.f. Fig.[Ij 6 ). The zero modes responsible for this 
periodicity are precisely the parafermionic zero modes 
Ti, Ai of order s = 3. 

We will show in Sec.|V]that the term proportional to k 
is a relevant perturbation. As the system size increases, 
smaller and smaller values of k suffice to produce an ap¬ 
preciable change in the periodicity of the supercurrent. 


V. STABILITY AGAINST LOCAL 
PERTURBATIONS AT ZERO TEMPERATURE 

As we saw in the previous section, gauge fluctuations 
change the topological edge structure of the parafermion 
chain. What is the stability of the modified edge against 
a generic perturbation? In order to gain some insight 
into this problem, in this section we determine numeri¬ 
cally the phase diagram of the gauged parafermion chain 
in the coupling space defined by h and k. Recall that 
the perturbation of strength h favors a trivial pairing of 
parafermions. For numerical investigations it is conve¬ 
nient to explicitly remove the gauge redundancy of the 
gauged parafermion chain. Hence in this section we work 




0 TT 27r Stt 47r Stt Ott 
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FIG. 1. Lowest energy levels of the Hamiltonian (I62II modeling a Ze parafermionic ring chain as a function of the enclosed flux 
The Ze parafermions are coupled to Ze gauge fields whose dynamics is restricted to Z 2 . (a) For static gauge fields with 
K = 0, the spectrum shows the 127r flux periodicity expected for Ze parafermions. ( 6 ) For dynamic gauge fields with k , = 0.2 J, 
the flux periodicity is reduced from 127r to a flux periodicity of Orr consistent with effective Z 3 parafermions. 


exclusively with the dual Hamiltonian of Eq. (ISTI) 

in the gauge invariant sector. 

Since energy levels, and quantum phase diagrams in 
particular, are insensitive to the distinction between 
parafermion and clock chains, let us emphasize once more 
the conceptual issue at stake. The gauged parafermion 
chain and the clock Hamiltonian of Eq. (HU), or its gauge- 
fixed version, are the same Hamiltonian from an spectral 
point of view. However, if one were to realize iLgp on a 
physical platform, then the representation Eq. (l3^ would 
be naturally suited for a fermionic platfor m^^d^ , and that 
of Eq. (l42|) a bosonic one, e.g., bosonic cold atoms. How¬ 
ever, for the Hamiltonian iLgp written as in Eq. (l42ll . it is 
easy to break some or all of its global symmetries explic¬ 
itly with perturbations of the form —hs H.c.). 

No gauge field is needed to effectively generate such 
terms. Hence, in a cold atom realization, the ground 
degeneracy of iLgp is fragilei^ In contrast, such perturba¬ 
tions are non-local in terms of parafermions and strictly 
absent in fermionic realizations, except, effectively, as 
gauge fields! 

In the following, we will focus on parafermions of or¬ 
der p = 6. Then it is possible to restrict the gauge field 
dynamics to with r = 2, 3 by taking a = 3, 2 respec¬ 
tively. We call the corresponding models the Z( 6 , 2) and 
Z( 6 , 3) chains. 


A. Critical exponents for Z(6,2) 

In the previous section, we have discussed that increas¬ 
ing h should eventually drive the system from a topolog¬ 
ical phase into a trivial phase. If we indeed have Z 3 
parafermionic zero modes for finite k, as suggested by 
the supercurrent numerical calculation of Sec. El this 
must be reflected in the universality class of the transi¬ 
tion. To confirm the presence of Z 3 parafermionic zero 
modes, we calculate the critical exponents of the transi¬ 


tion by means of a detailed finite-size scaling analysis. 

For the analysis of the system, we use the Hamilto¬ 
nian m in the physical sector = 0 which decom¬ 
poses into subspaces Hq to eigenvalue of the con¬ 

served quantity Vi. We assume that a single relevant 
operator associated with h will drive a transition from 
a symmetry-broken Z 3 phase into a trivial phase. Con¬ 
sequently, a scaling analysis along the lines of the Ising 
transition in absence of symmetry-breaking fields (in the 
bulk) applies 1 ^ While the bulk terms of Hg are invari¬ 
ant under the Z 3 symmetry, the operator vn and the 
operator + H.c) arising from the projection 

on Hg explicitly breaks this symmetry at the boundary. 
Although negligible in the thermodynamic limit, those 
symmetry-breaking boundary operators are expected to 
give rise to irrelevant scaling fields leading to increased 
finite-size effectsi^ 

For the numerical determination of the location of the 
critical point and the critical exponents, we use the den¬ 
sity renormalization group (DMRG) algorithm and vari- 
ationally determine the eigenstates of the Hamiltonian 
Eq. (1511) in the ground state sector <7 = 0, requiring for 
each variational state |' 0 ) with (approximate) eigenvalue 
A a precision of at least /\ < 10“^. 

We analyze systems with system size N = 16 , 32, 64 , 128 
sites and maximal internal bond dimensions D = 128. 

Let us start with a detailed finite-size scaling analy¬ 
sis along the line h = k. The similarity to the Ising 
case suggests exploiting the Binder cumulant Ul = 
1 — {ra ^)of the “magnetization per site” m = 
iyi+v\)/2N for a precise determination of the loca¬ 
tion of the critical point. Considering one irrelevant scal¬ 
ing field with renormalization group eigenvalue —uj < 0 
and expanding the scaling fields to lowest non-vanishing 
order, the Binder cumulant satisfies the scaling relation^ 

UL = ^{\h-K\N^/\u^N-^). ( 64 ) 

Here $ is a universal scaling function, v is the critical 
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hjJ hjj 


FIG. 2. Behavior of the Binder cumulants Ul = 1 — 
(m‘^)/3(m^)^, defined in terms of the magnetization m = 
for systems of different sizes N along the 
line h = K, in vicinity of the crossing points. 


FIG. 3. Behavior of the susceptibility x along the line h = k, 
for different system sizes N. Data points are not indicated 
for better visual clarity. 


exponent associated with the divergence of the correla¬ 
tion length, he is the critical point, and uq is the value 
of the irrelevant scaling held at the critical point. As be¬ 
fore, N denotes the system size. By expanding $ around 
the critical point (0, uo), we hnd that the crossing points 
hNi,N 2 of tho Binder cumulants of different systems with 
size Ni, N 2 scale as 


\hNi,N2 


hc\ (X 




1 - (A^2/iVi)~" 
1 - {N2/Ni)-y- 


(65) 


and thus converge to the critical point he as Ni goes to 
inhnity for hxed ratios Ni/N 2 . 

Figure Fig. [2] shows the behavior of the Binder cu¬ 
mulants for different system sizes N = 16, 32,64,128 and 
D = 128. We obtain the location of the critical point and 
the exponent —1/v — w using a Shanks transformation 
for the data tuples {Ni,N 2 ) = (16,32), (32,64), (64,128). 
Since we only have 3 data points available, we estimate 
the errors by varying he and observing when the data 
points ~ ^c| in a double-logarithmic plot devi¬ 

ate visibly from a line, yielding 1/zz -|- w Ri 1.69 and 
he = 0.38426(5). In absence of irrelevant scaling helds, 
the scaling form Eq. (IMl) predicts that the values of t/i 
of different system sizes N collapse onto a single curve 
as a function of z = (h — he)N^/‘^. The failure of a data 
collapse assuming the absence of irrelevant scaling helds, 
using 1/jz Ri 1.69, conhrms the non-negligible role of ir¬ 
relevant scaling helds for the scaling analysis. 

For the determination of the critical exponents ly, 
we consider the susceptibility 


y = — (to)^). (66) 

In absence of irrelevant scaling helds, y satishes the scal¬ 
ing relation 


(67) 



{h-he)N^^''/J 

FIG. 4. Data collapse of the susceptibility data shown in 
Fig. [ 3 ] with V = 0.832, y/iz = 1.731 determined from a 
minimization of the residuals Eq. dSSI in an interval of size 
Az = ±0.4 around the susceptibility maximum. The criti¬ 
cal point he is fixed at the value obtained using Binder cu¬ 
mulants. Away from the maximum, no complete data col¬ 
lapse is obtained. To the right of the maximum, the curves 
for different N approach the maximum Amax = 128 curve 
monotonously from below, whereas to the left of the max¬ 
imum, they approach the Amax curve monotonously from 
above. This systematic structure of the deviations suggests 
contributions due to irrelevant scaling fields, possibly associ¬ 
ated with the symmetry-breaking boundary operators in the 
Hamiltonian Hq, as a possible explanation. Data points are 
not indicated for better visual clarity, the data point spacing 
is roughly Az = 0.1 for all curves. 


where <1> is a universal scaling function and y is the critical 
exponent associated with the divergence of the suscepti¬ 
bility at the critical point. The behavior of the suscepti¬ 
bility across the transition is shown in Fig. [3] for different 
system sizes. According to the scaling form Eq. (1571) . 
the susceptibility values yA”'*''^'^ obtained for different 
system sizes collapse onto a single curve as a function 
of z = {h — hc)A^/^. Then, for given values of z and 
he obtained using Binder cumulants, we extract the crit- 


X = N^/'^^\h-he\N^/''). 
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{h-K)N^/''/J 

FIG. 5. Data collapse of the Binder cumulants Ul obtained 
using he, V determined from Binder cumulants and suscepti¬ 
bility. Data points are not indicated for better visual clarity, 
the data point spacing is roughly A 2 : = 0.1 for all curves. 
Insets show close ups of the curves in the marked regions, 
showing that the distance to the A^max = 128 curve system¬ 
atically decreases as N is increased. 



FIG. 6. Data collapse of the susceptibility curves obtained for 
the model along the line k = 0 with p = 2 which corre¬ 
sponds to an Ising model with symmetry-breaking boundary 
conditions, using the exactly known critical point hejJ = 1 
and the critical exponents = 1, 7 = 7/4. The collapse shows 
a structure similar to the collapse of the susceptibility data 
obtained along the line k — h for p = 6 that was shown in 
Fig. |4l with deviations from perfect data collapse away from 
the susceptibility peak. 


ical exponents v, ^jv from the susceptibility values by 
minimizing the quadratic differences 


z i<j 


( 68 ) 


between the rescaled susceptibility curves of sys¬ 

tems of size Ni- Performing the ht in a range of 
^zj J = ±0.4 around the maximum of the scaled suscep¬ 
tibility, which is the most prominent feature, we obtain 


the data collapse shown in Fig. |4]with ^jv = 1.731(2), 
V = 0.832(3). 

The error estimates are obtained by checking when 
data collapses obtained using values i/±Ai/, 7 /i/±A( 7 /z/) 
for the critical exponents, with v, 7 / 1 / as quoted above 
and Ai/, A{'y/v) some perturbation, show appreciable 
deviations from the data collapse of Fig. 0] While by 
choice of the ht region, the collapse shown in Fig. 2] is 
very good around the maximum, it shows systematic de¬ 
viations away from it. We note in particular that the 
distance to the maximal size iVniax = 128 curve is de¬ 
creasing monotonously as N is increased. Performing the 
Ht around other regions of the susceptibility curve yields 
data collapses with deviations that do not show a similar 
internal consistency. A data collapse of the Binder cu- 
mulant Ql obtained with he and v as determined above 
shows reasonable data collapse with deviations showing 
the same systematic behavior, see Fig. [5] 

The systematic behavior of the deviations suggests 
irrelevant scaling Helds due to the symmetry-breaking 
boundary operators in the Hamiltonian Hq as a plausible 
explanation. In order to show that this explanation is in¬ 
deed consistent, we have simulated the model m with 
p = 2 along the line k = 0, where it is just an Ising model 
with symmetry-breaking boundary conditions. The col¬ 
lapse of the susceptibility curves of systems with different 
sizes Ni, using the exactly known values for the critical 
point he/ J = 1 and the critical exponents ^= 1,7 = 7/4, 
is shown in Fig. [51 The data shows good collapse around 
the susceptibility maximum and deviations away from 
it, i.e., a structure similar to the collapse discussed pre¬ 
viously. This shows that symmetry-breaking boundary 
operators may indeed be the cause of the irrelevant scal¬ 
ing Helds and strongly supports the assumption that our 
choice of Ht region around the susceptibility maximum al¬ 
lows us to extract the critical exponents with a negligible 
systematic error. 

Summing up, we Hnd v = 0.832(3), ^/v = 1.731(2) 
for the critical exponents of the transition along the line 
h = K. This result is in good agreement with the ex¬ 
act values v = 5/6 = 0.83, 7 /z/ = 26/15 = 1.73 ex¬ 
pected for the three-state Potts modeli^ We thus nu¬ 
merically conHrm the transition from a Z 3 to a trivial 
phase that was suggested by the structure of the dual 
Hamiltonian (1511) and the periodicity of the supercurrent 
discussed in Sec. lYl 


B. Phase diagram 

According to our Hndings for the Z( 6 , 2) model, the 
crossings of the Binder cumulants indicate the location 
of the critical point up to small Hnite-size corrections. In 
fact, this is expected to hold regardless of whether the I 12 
subgroup or the Z 3 subgroup of the gauge Helds is dy¬ 
namic and we will also use Binder cumulants to study the 
phase boundaries for the !Z(6,3) model. In this model, 
the parafermion edge modes of degree p = 6 are trans- 
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FIG. 7. Location of the critical line for the transition from the 
nontrivial phase (to the left of the line) to the trivial phase 
(to the right of the line) for the Z(6,2) and Z( 6 , 3) model. 
The critical line is determined from the crossing points of 
the Binder cumulants Ul for system sizes (A^i, N 2 ) = (32, 64) 
and the Binder cumulants are calculated numerically using 
the DMRG algorithm with internal bond dimension D — 64. 
Data points are marked as small crosses and linearly interpo¬ 
lated for visual guidance. The large cross indicates the posi¬ 
tion of the point where we have determined the critical expo¬ 
nents u = 0.832, — 1.731 consistent with a Z 3 phased 

by means of a more detailed hnite-size scaling analysis. The 
inset shows a closeup of the region around k = 0 , showing 
that in both cases, the critical lines of the models end at the 
same point at k = 0 . 


muted into Majorana edge modes. 

To determine the phase boundaries of both models, 
we use the crossings of the Binder cumulants of system 
sizes {Ni,N 2 ) = (32,64) with internal bond dimension 
D = 64. For the !Z( 6 , 2) model, we investigated the cross¬ 
ing along lines of constant k, for k/J > 0.2, whereas for 
k/J < 0.2, we focused on lines of constant h in order 
to to ensure that the critical line is hit at approximately 
right angle. We proceeded similarly for the Z( 6 , 3) model. 
The resulting phase diagrams of the Z( 6 , 3) and Z( 6 , 3) 
gauged parafermion chains are shown in Fig. [T] There is 
a topologically non-trivial Z 3 phase for the Z( 6 , 2 ) model 
and a topologically non-trivial Z 2 (Majorana) phase for 
the Z( 6 ,3) model in an extended region of the phase dia¬ 
gram. Notice that the range of Z 3 criticality in coupling 
space encompasses a line rather than a point as for the 
standard Z 3 parafermion chain. This observation may 
prove useful in connection to a recent blueprint for ob¬ 
taining Fibonacci anyons;— we will come back to this 
point in the outlook. 

Finally, let us emphasize that our numerical findings 
indicate that the two critical lines share a common end¬ 
point marking the end of the critical phase of the Zg 
parafermion chain. Obviously, our results are not to¬ 
tally conclusive due to the limited system sizes and the 
lack of detailed finite-size extrapolations in that area 
of the phase diagram. However, we can offer a quali¬ 
tative explanation of this numerical observation. It is 


a well-investigated fact that the perturbed, ungauged 
parafermion chain (a.k.a, the clock model) displays a 
critical phase for p > 5, see for example Ref. for a 
recent, comprehensive discussion. Imagine now for the 
sake of the argument that the transition line for the 
Z( 6 , 2) model were to meet the critical phase of the Zg 
parafermion chain somewhere inside that phase. Then, 
from that point on, that is, for larger values of h, the 
chain would effectively behave as a Z 2 chain, since the 
only edge inodes left would be the critical Majorana 
modes. But a Z 2 chain does not support a critical phase. 
Hence it must be that the critical lines meet at the end 
of the critical phase. Incidentally, this determines com¬ 
pletely the boundary of the critical phase, since the start¬ 
ing point is a function of the final point by self-duality. 


VI. SUMMARY AND OUTLOOK 

Topological zero-energy boundary modes are typically 
stable against generic perturbations, and it is precisely 
this feature that makes them attractive for quantum in¬ 
formation processing. However, it is also this feature 
that makes them hard to mold in a controlled fashion. 
Extending our previous work on the effect of phase slips 
in the Majorana chain of Kitaev, we have argued in this 
paper that for topological boundary modes associated to 
a global protecting symmetry, it is possible to modify the 
topological edge structure in a controlled fashion by en¬ 
gineering a gauge field into the system. The quantum 
fluctuations of the gauge field act as a relevant perturba¬ 
tion without symmetry breaking, and permit to split the 
ground degeneracy of the system, partially or completely, 
depending on the designed properties of the gauge fluctu¬ 
ations. If the gauge field is engineered, by allowing only 
restricted gauge fluctuations, so that the ground degen¬ 
eracy is only partially split, then the topological edge 
modes are modified accordingly, in a controlled and pre¬ 
dictable fashion. There is absolutely no risk of driving the 
system to a trivial phase purely by restricted gauge fluc¬ 
tuations, unless the gauge field is explicitly designed to 
do so by allowing for unrestricted gauge fluctuations. It 
would be interesting to recast these results in terms of the 
group cohomology classification of symmetry protected 
topological quantum orders and bosonic anomalies.— 

We have illustrated our ideas with gauged parafermion 
chains. For generic chains, we performed a symmetry and 
duality analysis conhrming our theoretical predictions. 
To obtain a more refined picture of the the transmutation 
of the edge modes and the correlation with critical phe¬ 
nomena, we investigated numerically the Zg parafermion 
chain coupled to a Z 2 -like or a Zg-like gauge held. That 
is, the dynamics of the physically natural Zg gauge held 
was restricted to achieve a Z 2 -like or Z 3 -like effect. Ac¬ 
cording to our general picture, the Z 2 -like held should 
transmute the Zg edge modes of the parafermion chain 
into Z 3 edge modes, and we showed numerically that this 
is the case by computing the period of the supercurrent in 
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a parafermion ring junction as a function of gauge fluctu¬ 
ations. We also obtained the phase diagram of the gauged 
parafermion chain as a function of gauge fluctuations and 
a perturbation driving the system to a trivial phase. In 
this two dimensional phase diagram, the transition be¬ 
tween the gauge-driven phase with Z 3 edge modes and 
the topologically trivial phase occurs on a critical line in 
the universality class of the Z 3 parafermion chain. From 
this point of view, the Zi 2 -gauged Ze parafermion chain is 
indistinguishable from the ungauged Z 3 chain. The gen¬ 
eral picture is the same for the Z 3 -gauged Zg parafermion 
chain. In this case, the edge modes are Majoranas. 

Let us mention in closing two potential practical ap¬ 
plications of our work. One possibility is to exploit our 
ideas to create topological qutrits out of Zg parafermions. 
Qutrits have practical advantages over qubits, not be¬ 
cause they can provide qualitatively faster algorithms, 
but because they can polynomialy shorten the length of 
an algorithm in the circuit model of quantum computa¬ 
tion. Of course, the same holds for six-state quantum 
bits, but here we have face the problem that very little is 
known about quantum software design with many-leveled 
logic. Topological qutrits strike a nice balance from this 
point of view. 

Another possibility is to use a gauged Zg parafermion 
chain as a replacement for a Z 3 chain. The Z 3 


parafermion chains is one of the key building blocks in 
a recent blueprinti^ for realizing Fibonacci anyons. One 
difficulty of that blueprint is that, as a building block 
in this setup, all the chains in a two-dimensional stack 
must be tuned to their critical point. By contrast, our 
realization of the Z 3 chain out of the Zg chain is critical 
on a line rather than a point, and the whole line is in 
the universality class of the Z 3 chain. Moreover, the Zg 
chain is in principle realizable out of the (doubled) 1/3 
fractional quantum Hall state, which is the most stable 
fractional plateau. We expect that the Z 2 gauge field for 
this application might also be engineered out of super¬ 
conducting phase slips, but the situation is not as simple 
as for the Majorana chain. The mesoscopic details of en¬ 
gineering the required gauge fields will be the subject of 
an upcoming publication. 
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